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ABSTRACT 


This  thesis  presents  a  discussion  of  the  problems 
involved  in  estimation  of  Circular  Error  Probable  (CEP) . 
Several  estimators  are  compared  through  simulation,  under 
two  models;  the  power  model  and  the  Rayleigh  model.  Several 
measures  of  effectiveness  are  calculated  for  each  of  the 
competing  estimators.  It  is  found  that  maximum  likelihood 
estimation  based  on  the  power  distribution  performs  well 
for  "heavy  tailed"  distributions;  the  Rayleigh  unbiased 
estimator  performs  well  for  most  other  situations. 
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INTRODUCTION 


I  . 

Evaluation  of  a  weapons  system  is  frequently  an  expen¬ 
sive  and  time  consuming  task.  Only  relatively  few  tests 
may  be  allowed  because  of  the  high  cost  of  testing.  It  is 
possible  that  a  good  weapon  system  could  be  rejected  because 
of  inefficient  utilization  of  the  small  amount  of  data 
available.  Hence,  efficient  use  of  the  limited  data  is 
required  to  reduce  this  risk. 

Several  authors  have  considered  the  measurement  of 
weapon  delivery  accuracy  using  Circular  Error  Probable 
(CEP) .  This  measure  of  accuracy  is  defined  to  be  the  radius 
of  a  circle  centered  at  the  target  which  on  the  average  con¬ 
tains  fifty  percent  of  impact  points,  that  is,  the  median 
radial  miss  distance.  A  functional  relationship  between 
CEP  and  the  distribution  parameters  is  given  by  the 
following  integral: 

CEP 

F(CEP)  =  /  f(r)dr  =  0.5  , 

0 

where  f(r)  is  the  probability  density  of  radial  miss 
distance . 

Various  parametric  models  for  the  distribution  of  impact 
points  about  the  target  have  been  used  in  connection  with 
such  estimations.  A  discussion  of  these  models  is  given  by 
Jordan  [1] .  But,  deficiencies  in  the  models  have  been 
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noted.  Specifically,  they  tend  to  provide  poor  fit  to  the 
"upper  tail"  of  the  radial  distributions.  That  is,  they 
provide  less  than  satisfactory  explanation  of  impacts  ob¬ 
served  far  from  the  target.  A  commonly  used  model  is  the 
"Rayleigh  model",  in  which  impacts  in  the  target  plane  are 
assumed  to  be  distributed  as  the  unbiased  circular  normal 
distribution.  One  might  anticipate  that  Rayleigh-based 
CEP  estimators  might  not  perform  adequately  with  a  large 
proportion  of  actual  test  data.  Several  potential  competitors 
to  Rayleigh  CEP  estimators  were  therefore  developed  for 
comparison . 

Barr  and  Jordan  [2]  developed  and  proposed  the  power 
distribution  model  which  is  based  on  an  infinite  mixture  of 
Rayleigh  radial  miss  distance  distributions.  This  distri¬ 
bution  has  properties  compatible  with  the  interpretation  of 
a  in  the  underlying  Rayleigh  distribution  model;  they  sug¬ 
gested  the  use  of  maximum  likelihood  estimates  of  the  param¬ 
eters  of  the  power  distribution  which  in  turn  provide  a 
maximum  likelihood  estimate  (MLE)  of  CEP. 

The  jackknife  is  a  technique  which  v;as  proposed  by 
Quenouille  [10]  for  reducing  bias  in  estimation  problems. 

We  tried  the  jackknife  procedure  to  estimate  CEP  under  the 
Rayleigh  model. 

The  sample  median  was  also  used  to  estimate  CEP.  Since 
CEP  is  defined  as  the  median  radial  miss  distance,  a  natural 
estimator  for  CEP  is  the  sample  median. 
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To  measure  the  performance  of  the  estimators,  a  variety 
of  data  sets  from  known  distributions,  such  as  uniform, 
Rayleigh,  pov;er  and  Weibull  were  generated  by  simulation. 
Samples  from  the  known  distributions  were  generated  to 
test  the  capability  of  the  estimators  to  handle  widely 
disparate  situations. 
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II.  MAXIMUM  LIKELIHOOD  ESTIMATION  BASED 
ON  THE  POVJER  DISTRIBUTION 

Barr  and  Jordan  [1],  [2]  have  given  a  discussion  of 

CEP  accuracy  measurement  under  a  pov;er  distribution  model. 

They  showed  that  estimating  the  parameters  of  the  power 

distribution  by  the  method  of  moments  may  not  be  adequate; 

moreover  it  is  possible  only  when  the  ratio  of  the  squared 

—2  2 

sample  mean  and  variance  r  /s^  does  not  exceed  7t/4  -  tt. 
Consequently,  they  suggested  use  of  the  maximum  likelihood 
(M.L.)  method. 

Analytically  closed  expressions  for  the  M.L.  estimators 
cannot  be  obtained,  the  solution  of  normal  equations  yielding 
the  parameter  estimates  requires  numerical  techniques.  In 
this  section  we  examine  a  numerical  method  to  obtain  the 
M.L.  estimates  of  the  parameters  of  the  power  distribution, 
which  in  turn  provide  a  M.L.  estimate  of  CEP.  These  estimation 
procedures  require  use  of  a  computer. 

A.  MATHEMATICAL  FORMULATION 

The  probability  density  function  of  the  radial  miss 
distance  R  having  a  power  distribution  with  parameters  k  and 
?  is 


fj^(r  ;k,  ;) 


2<r 


(- 

c 


+  r 


K  +  1 


r  >  0 ,  K ,  ^  >  0 
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A  C.D.F.  for  radial  miss  distance  is 


Fj^(r ;  K, 


1  -  (- 


+ 


< 


r>0,  k:,c>0  (1) 


The  CEP  is  given  by 

1  1 

CEP  =  =  [(2^  -  1)^]^ 


Because  radial  miss  distance  is  alv/ays  non-negative, 
one  can  equivalently  consider  squared  radial  miss  distance. 

In  the  case  of  the  power  model,  this  leads  to  some  notational 

.  .  .  .  2 
simplification.  Let  Y  denote  the  squared  miss  distance  R  . 

Then  from  (1) , 

F^Cy)  =  1  -  (1  +  ;  y  >  0  ,  K,^  >  0 

(2) 

fY(y)  =  f(l  +  ?  y  1  0  ,  >  0 

Here  ^  and  k  are  parameters  to  be  estimated  from  observed 
impact  data.  The  relationship  between  CEP  and  these 
parameters  is 

1 

CEP^  =  Fy~^(|)  =  [(2^^  -  1)^] 


The  method  of  maximum  likelihood  is  based  upon  the 
likelihood  function  ({)(y)  which  is  defined  as  the  joint 
density  function  of  a  sample  of  size  n;  i.e.. 
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n  n  n  Y  •  _  t 

4>(y)  =  TT  fyCyJ  =  %  TT  (1  + 

i=l  c  i=l  ^ 

since  it  is  easier  to  deal  with  sums  rather  than  products 
when  maximizing,  <()(y)  is  transformed  into 

L(y)  =  log  (J)(y) 

n 

=  n  log  K  +  nxlog  ^  -  (k+1)  Z  log(^+y.) 

i=l  ^ 

There  is  no  loss  of  generality  here  because  the  maximum  of 
a  positive  function  occurs  at  the  same  point  as  the  maximum 
of  the  logarithm  of  the  function.  The  normal  equations  are 
then 


m(K , ;) 


_  3L(y) 
8k 


n 

^  +  nlog  C -  Z  log(^+y.) 

i=l  ^ 


0 


and 


g(Kf  C) 


=  9L(y) 
9C 


nK 


(K  +  1) 


n 

Z 

i=l 


1 


^  +  Y; 


0 


The  solution  of  m(S,t)  =0  and  g(K,c)  =0  gives  the 
parameter  M.L.  estimators  k,  As  mentioned  before,  in 

the  present  case,  this  requires  numerical  techniques. 


B.  NUMERICAL  SOLUTIONS:  NEWTON'S  METHOD 

This  method  is  a  Newton's  iteration  technique  and  the 
concept  and  procedure  are  discussed  in  [3].  Newton's  method 
provides  a  statistical  criterion  for  stopping  the  iterations. 
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The  pure  form  of  Newton's  method  is  given  by  the 


following  iteration  scheme: 


/\  ^  T  m 

■  tf’(Xj)  ]  Vf  (Xj)  ^ 


where 


_  /S  ^  -m 


f(x.)  =  [ii45±,iL(Xl, 


3k 


3? 


First  we  compute 


Let 


9L(y)  ^ 

/s 

3< 


n 


n 


-  +  n  log  ^  -  E  log  (^+y.) 
K  i=l  ^ 


3L(y)  ^ 

/N 

3C 


UK 


n 


-  (k+1)  E 


i=i  c  +  y. 


3^L(y) 

9k 


-n 

-2 

K 


3^L(y)  ^  n  _  ^  1 

3k3c  ;  i=l  ^  +  y^ 


3^L(y) 


- 


n 


+  (k+1)  E 


i=i  (C  +  y.) 


D  =  3^L(y)  .  3^L(y)  _  ( i^MxL)  2 

3k  3c  3k3c 
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If  D  were  identically  equal  to  zero,  this  method  would  not 
be  applicable.  Now  take 


-1  '' 

F  (X)  = 


9^L(y)  _  8^L(y) 

3c  3k3c; 

9^L(y)  3^L(y) 

9k:9c 


D 


In  order  to  develop  the  termination  criterion,  consider 
the  statistic  for  (for  reference  see  [4]) 

B  =  Vf (X.)f“^(X.) Vf (X.)”^ 

1  1  1 

B  is  asymptotically  a  Chi-square  r.  v.  with  2  degrees  of 

freedom,  therefore  it  might  be  reasonable  to  continue  the 

iteration  until  the  statistic  B  becomes  less  than  the  100a 

percentile  point  of  the  chi-square  distribution  with  two 

2 

degrees  of  freedom,  x  • 


14 


III.  MAXIMUM  LIKELIHOOD  ESTIMATION  BASED 
ON  THE  RAYELIGH  DISTRIBUTION  MODEL 


Suppose  impacts,  measured  in  a  suitable  plane,  have 

range  (X)  and  deflection  (Y)  components  of  impact  that  are 

jointly  uncorrected  bivariate  normal  distributed  with  mean 

2 

at  the  target  and  common  unknown  variance  a  in  both  direc¬ 
tions.  That  is,  for  the  present  development,  it  will  be 
assxamed  that  the  joint  density  of  miss  distances  about  the 
target  is  given  by 

4:1  \  2,-1  -(x^+y^)/2o^ 

f(x,y)  =  (2tt0  )  e'  ^  '  ,  -a><x,y<“. 

It  is  easy  to  derive  the  density  of  radial  miss  distance. 
1/2 

Let  R  =  (X+Y)  '  and  transform  to  polar  coordinates  to 
obtain 

T  2Tr  r  - 2 

F  (r)  =  P[R£r]  =  -- ---y-  /  /  R  e  dR  d0 

2Tra  0  0 

-  2 
r 

",  2 

=  1  -  e  °  ;  r  ^  0,  0  >  0  . 

This  distribution  of  the  radial  miss  distance  is  referred 
to  as  either  the  radial  normal  distribution  or  the  Rayleigh 
distribution . 

Then  the  density  function  of  the  radial  miss  distance, 
fj^(r)  is 
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f^(r) 


r  2g 

~7  ® 
a 


r  >  0,  o  >  0 


To  find  the  relationship  between  CEP  and  o,  we  calculate 


CEP  =  =  y2  In  2  a  :  1.1774  a  , 


which  is  a  well  known  relation.  Hence,  if  we  want  to  estimate 
CEP  with  this  model,  first  we  might  estimate  the  population 
standard  deviation.  The  maximum  likelihood  estimator  for 
o  is  easily  found  to  be 


a 


2  2 

i=l  ^ 


f 


where  r^,  ...,  r^  is  the  sample  of  observed  radial  miss 

distances.  This  estimator  has  a  slight  bias,  however.  An 
adjusted  unbiased  estimator  for  CEP  is 


CEP  =  1.1774  (Jn  -lilli - )  g,  (3) 

r{|(2n+l)) 

where  the  quantity  in  parenthesis  is  an  unbiasing  factor 
for  a.  It  has  been  shown  by  Chapman  and  Robbins  [6]  that 
CEP  =  1.1848  a  for  sample  size  20  is  a  best  estimator  of  CEP 
under  the  normality  and  other  assumptions  stated  above. 
However  other  estimators,  while  not  as  efficient  as  the  best 
estimator,  may  have  more  robustness  under  variations  in 
the  assumed  distribution  of  radial  miss  distance. 
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The  term  "jackknife"  has  been  used  by  Miller  [9]  and 

Quenouille  [10]  to  describe  a  method  of  modifying  a  biased 

estimator  to  bring  it  closer  to  an  unbiased  condition. 

From  (3)  we  see  the  M.L.  CEP  estimator  has  a  bias  term  of 

order  — .  We  are  thus  led  to  try  to  examine  the  jackknife 

/n 

method  to  reduce  this  bias  in  the  estimator.  Brillinger [7 ] 
has  provided  a  discussion  of  the  asymptotic  properties  of 
estimators  obtained  with  the  jackknife  method  when  it  is 
applied  to  M.L.E's. 

The  jackknife  procedure  is  based  on  dividing  the  data 
into  groups,  obtaining  estimates  from  combinations  of  the 
groups,  and  then  averaging  these  estimates.  We  apply  this 
to  our  problem,  as  follows.  The  maximum  likelihood  estimator 
for  o  IS 


-2 

o 


r  . 
D 


2 


Divide  these  data  into  m  groups,  where  r^,  ...,  r^  is  the 

sample  of  n  observed  radial  miss  distances  and  the  size  of 
each  subgroup  is  k.  First  estimate  the  parameter  o  based 
on  all  data 


2 


1 

^  E  r  . ' 
2n  3 


'^2  th 

next  estimate  o  based  on  all  data  except  those  in  the  i 

group 
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I 


-  2  1^2  ... 

“-i  =  ^  ^  =  1,  2,  . 


•  t  m  • 


Then  compute  the  pseudo  values 


^2  ''2  ''2 

o.  =  ma  -  (m-l)a  .  ;  i  =  1,  2,  ..., 

1  O  - 1 


Finally,  compute  the  average  of  the  pseudo  values 

^  o 

I  O. 


m 


and  estimate  the  CEP  by 


CEP 


^2 

In  2  a 


m. 
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IV.  THE  SAMPLE  MEDIAN 


Since  CEP  is  defined  as  the  median  radial  miss  distance, 
a  natural  estimator  for  CEP  is  the  sample  median  (m) .  By 
definition;  the  median  is  the  50^  percentile  of  the  distri¬ 
bution,  i .  e .  , 

P(R  £  m)  ^0.50  and  P(R  ^  m)  ^  .50 

Since  F(CEP)  =  0.50  and  R  is  assumed  to  have  a  continuous 
distribution. 


-11 

CEP  =  m  =  F  (i) 

Suppose  that  the  n  radial  miss  distance  observations  of 
the  sample  are  arranged  in  order  of  magnitude,  ^(2)' 

...,  shall  consider  only  the  estimates  of  the 

population  median  defined  as  follows: 

(a)  If  the  sample  size  is  odd,  and  n  =  2p+l,  the  median 
(m)  is  taken  as  the  (p+1) ^  value. 

(b)  If  the  sample  size  is  even,  and  n  =  2p,  the  median 
(m)  is  taken  as  the  midpoint  between  the  p  and 
(p+1)^^  values  or 


1 

7 


(r 


(P) 


+ 


"^(P+I) 


)  . 
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The  distribution  of  the  median  will  tend  to  be  symmetri¬ 
cal,  and  should  be  very  nearly  normal  in  almost  all  popula¬ 
tions.  The  distribution  of  the  median  in  the  case  of  small 
samples  from  a  normal  population,  whether  the  sample  size 
is  odd  or  even,  will  tend  to  be  normal  [Hoji  [11]) .  The 
normal  form  is  approached  very  rapidly.  The  value  of  the 
standard  deviations  of  the  median  for  n  even  and  n  odd 
approach  the  same  limit. 


where  a  is  the  standard  deviation  in  the  population  sampled. 

We  use  the  order  statistics  to  examine  the  sample  median 
under  the  Rayleigh  distribution  model  (see  Inselmann  [12]), 
and  find  its  bias  and  variance  so  as  to  enable  one  to  ma)ce 
a  comparison  with  other  estimates.  If  n  is  even,  the  dis¬ 
tribution  of  sample  median  is  given  by  the  following: 


g(r(0  5n))  =  {n!/{n-[0.5n])  !  (  [0.5n]-l)  !}f  {r)F{r  JQ 


[0. 5n]-l 


j  ^n- [0. 5n] 


For  convenience  we  let  a  =  0.5n,  and 


C  =  n!/(n-[0.5n])  !  (  [0.5n]-l)  !  =  (2a)  './a !  (a-1)  ! 
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Then  we  may  write  the  distri- 


so  the  sample  median  is  r^. 
bution  of  sample  median  under  the  Rayleigh  distribution  model, 


2  2  _  2 
g(ra)  =  c[-^exp{ — ]  [1-exp  (-- ^-2-)  ]  ^  [ exp  ( — ® 

o  2o  2a  2a 

2  2 
r  r  I  ^  (s+1)  \  1  n  !  ^  \ia~l 

=  c  (exp(-  2 — )][1  -  exp{ - j)  ] 

a  2a  2a 


We  will  now  consider  computation  of  the  mean  of  ra. 


E(ra)  =  c  /  ^(exp(-^— ]  [1-exp (-  5-)]^  dr 

0  a  2a  2a 


f”  ,  r^(a+l),i  ,  ,.K,a-l,  r  ,  r^.,Kj 

=  c  /  -^[exp(-  - - i-)  ]  Z  (-1)  {  )  [exp( - g’)  ^ 

0  a  2a  K=0  2a 


,  .,K,a-l,  “  r^  ,  r^(a+l+K),  . 

=  c  Z  (-1)  {  y  )  /  “T  exp( - i - 5- — -)  dr 

K=0  ^  0  a  2o 


and  make  the  following  transformation 
y  =  r^ (a+l+K) /2a^  , 
r  =  a  J  2y/ (a+l+K)  , 

dy  =  [r (a+l+K) /a^ ]  dr  . 
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To  obtain, 


E(ra)  =  co/T  Z  (-1)  ^  (\^)  /  (a+l+K)  ^  /  ./Y  e  ^  dy 

K=0  ^  0 


3  a-1  (-1)^(\^) 

=  car(4)-/T*  E  - 

K=0  (a+l+K)  / 


We  can  adjust  the  bias  term  and  using  the  same  transformation 
as  above  we  may  now  find  the  second  moment 


2 

„  “>  3  2,  a-1  ,  ,>K,a-l>  ,  ,  t  ,iK  , 

E(r  2)  =  c  /  E  ^  K  )texp(-^)]  dr 

^  0  2a^  K=0 

2  a-1  (-l)^(^^^) 

=  2co  Z  - ^ 

K=0  (a+l+K)^ 


The  variances  of  this  sample  median  can  now  be  computed 
using  the  well-known  formula 


Var  (r^)  =E(r^^)  -  [E(r^)]^  . 
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V.  SIMULATION  NJITHODOLOGY 


The  previous  sections  have  been  concerned  with  a 
development  of  different  estimators  for  the  median  of  the 
distribution  of  radial  miss  distances  obtained  under  a  set 
of  fixed  impact  conditions.  An  analytical  comparison  is 
difficult  to  perform  for  these  estimators,  so  computer 
simulation  was  used.  This  section  summarizes  the  different 
models  and  estimators  discussed  in  the  previous  sections, 
and  includes  an  analysis  of  the  results  obtained  from  simu¬ 
lations.  Although  the  sample  problems  do  not  represent  actual 
weapon  test  results,  an  attempt  has  been  made  to  generate 
data  from  a  wide  variety  of  "realistic"  distributions. 
Therefore  this  simulation  analysis  should  show  relationships 
between  our  estimators  of  CEP,  and  indicate  their  relative 
robustness,  and  hence  potential  for  use  in  a  variety  of 
applications.  The  alternative  estimators  are  as  follows: 

/V 

CEP  1:  maximum  likelihood  unbiased  Rayleigh  estimator 

CEP  2:  maximum  likelihood  Jackknife  Rayleigh  estimator, 
with  subgroup  sample  size  k  =  1 

CEP  3;  maximum  likelihood  power  estimator  under  Newton's 
method 

/N 

CEP  4:  sample  median  estimator 

For  the  purpose  of  our  simulation  comparison  of  these 
four  CEP  estimators,  we  generated  samples  of  radial  miss 
distances  (or,  in  some  cases,  squared  radial  miss  distance) 


23 


from  several  parent  distributions:  Rayleigh,  uniform,  power 
and  Neibull  distributions.  Throughout,  we  assumed  a  sample 
size  of  20  impact  observations.  The  method  of  generation 
was  by  C.D.F.  inversion  of  generated  uniform  (0,1)  values. 
For  each  parent  distribution  (hereafter  called  a  "case") 
we  replicated  estimation  of  CEP's  500  times.  Summaries  of 
the  sample  distributions  of  500  values  of  5  CEP's  were 
developed . 

Eight  cases  were  considered;  in  each,  parameter  values 
were  used  which  gave  a  true  CEP  of  100.  The  cases  are  as 
follows  (U  is  a  uniform  (0,1)  random  variable  throughout): 

Case  I.  (Rayleigh  distribution) 

2 

F  ^  (y)  =  1  -  e  with  =  10.^/21n2 

CEP  =  In  2  a 

Generated  value:  Y  =  -2a^  In  U  =  -14426.95  In  U 

Case  II.  (Uniform  distribution) 

FR(y)  =  y  ,  with  y  e  (0,200). 

CEP  =  100 

Generated  value  Y  =  200  U  . 
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Case  III.  (Power  distribution) 


Fj^^(y)  =  1  -  (1  +  |)"^ 


y  >  0 


>  0 


with  K  =  1  and  ^  =  10^ 

1  1 
CEP  =  (C  (2'^  -  1)  ]  ^ 


_1 

Generated  value  Y  =  (U  -  1)^  =  10^  -  1) 


Case  IV.  (Translated  Rayleigh  distribution) 

P  2(y)  =  1  -  ^-i^-l00)/2a^ 

K 

where  2a^  =  14282.6809 

1 

CEP  =  [100  +  2a^  In  2]^ 


2 

Generated  value  Y  =  -  2a  In  U  +  100 


Cases  V-VIII.  (Weibull  distributions) 

w 

V 

2 

Fr  (y)  =  1  -  e 
CEP  =  (v  In  2)^^^ 


Generated  value  Y  in  the  table  below 
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Four  sets  of  parameters  (w,v)  were  selected  to  give 
CEP  =  100  and  to  provide  distributions  with  varying  upper 
tail  thickness.  One  case  (with  w  =  1)  gives  an  exponential 
distribution  with  mean  v,  and  this  corresponds  to  the  Rayleigh 
case.  Case  I  above  may  be  considered  also  as  a  Weibull  case. 
The  parameter  pairs  used,  in  order  of  increasing  upper  tail 
size,  are  given  in  the  following  table. 


1 


Case 

w 

V 

Y  = 

(-vlnU) 

V 

2 

10®/ln2 

Y  = 

(-144260504.1  In  U) 

2 

VI 

1.5 

10^/ln2 

Y  = 

(-1442695.04  In  U) ^ 

I 

1.0 

10'^/ln2 

Y  = 

-4426.95  In  U 

VII 

.5 

10^/ln2 

Y  = 

(-144.27  In  U) ^ 

VIII 

.1 

10*'^/ln2 

Y  = 

(-3.623886  In  U) 
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VI.  SIMULATION  RESULTS 


As  mentioned  above,  we  tested  four  CEP  estimators  from 
eight  different  cases.  Major  comparison  measures  include 
bias,  variance  and  mean  square  error.  Others  are  included 
in  the  tables  that  follow. 
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CASE  I 


RAYLEIGH  DISTRIBUTION 


y  =  -  14426.95  In  U 


* 


MOE 

CEP  1* 

CEP  2 

CEP  3 

CEP  4 

MEAN 

99.93 

99.94 

95.96 

107 

STD  DEV 

11.66 

11.68 

12.22 

38.8 

MSE 

136.03 

136.37 

165.77 

1549 

GED  MEAN 

99.25 

99.25 

95.17 

100 

VARIANCE 

136.03 

137.37 

149.45 

1500 

COEFF  VAR 

0.12 

0.12 

0.13 

0.36 

SKEWNESS 

0.15 

0.14 

0.049 

0.35 

KURTOSIS 

0.10 

0.085 

0.307 

0.34 

RANGE 

71.38 

70.87 

83..60 

252 

MINIMUM 

67.25 

67.29 

53.74 

17.7 

.10  QUANTILE 

84.82 

84.80 

80.99 

58.4 

.25 

92.12 

92.09 

87.67 

82.5 

.50 

99.63 

99.63 

95.96 

105 

.75 

107.49 

107.34 

104.00 

132 

.90 

114.83 

114.83 

111.32 

156 

MAXIMUM 

138.63 

138.16 

137.34 

269 

CEP  1:  maximum  likelihood  unbiased  Rayleigh  estimator. 

CEP  2:  maximum  likelihood  Jackknife  Rayleigh  estimator, 
with  subgroup  sample  size  k  =  1. 

CEP  3;  maximum  likelihood  power  estimator  under 
Newton's  Second  method. 


CEP  4;  sample  median  estimator. 
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CASE  II 


UNIFORM  DISTRIBUTION 


y 

/N 

=  (200  U) 

/N 

MOE 

CEP  1 

CEP  2 

CEP  3 

CEP  4 

MEAN 

98.61 

96.52 

91.24 

102 

STD  DEV 

9.87 

9.73 

17.08 

41.7 

MSE 

109.32 

107.03 

368.63 

1744.0 

GED  MEAN 

96.09 

96.01 

* 

91.5 

VARIANCE 

97.83 

94.92 

291.89 

1740 

COEFF  VAR 

0.10 

0.10 

0.19 

0.40 

SKEV7NESS 

-  0.21 

-  0.20 

-  1.83 

0.005 

KURTOSIS 

0.61 

0.60 

5.54 

-  0.56 

RANGE 

67.28 

66.35 

125.60 

191 

MINIMUM 

59.52 

59.97 

0.0 

3.75 

.10  QUANTILE 

84.19 

84.25 

72.1 

45.7 

.  25 

90.10 

90.13 

85.30 

742 

.50 

96.84 

96.74 

94.26 

102 

.75 

102.96 

102.75 

101.52 

132 

.90 

108.66 

108.40 

107.62 

158 

MAXIMUM 

126.80 

126.32 

125.60 

194 

/s 

*  some  CEP  = 

0.0 
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CASE  III 


P0V7ER  DISTRIBUTION 


Y  =  •  (i  -  1) 


MOE 

CEP  1 

CEP  2 

CEP  3 

CEP  4 

MEAN 

209.61 

229.33 

100.71 

113 

STD  DEV 

190.15 

257.90 

21.07 

46.0 

MSE 

2460 

83000 

444.36 

2279 

GED  MEAN 

177.76 

185.62 

98.60 

103 

VARIANCE 

3615.71 

66515.7 

443.86 

2110 

COEF  VAR 

0.90 

1.12 

0.21 

0.40 

SKE^VNESS 

6.39 

7.32 

0.84 

0.76 

KURTOSIS 

57.54 

72.76 

2.05 

1.48 

RANGE 

2264.4 

3220.29 

158.52 

337 

MINIMUM 

74.59 

75.08 

48.32 

17.2 

.10  QUANTILE 

108.29 

108.79 

76.88 

57.1 

.  25 

126.07 

127.96 

86.19 

82.5 

.  50 

161.63 

166.20 

99.25 

107 

,75 

220.09 

228.03 

112.71 

140 

.90 

331.77 

364.51 

124.77 

172 

MAXIMUM 

2339.01 

3295.38 

206.84 

354 
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CASE  IV.  TRANSLATED  RAYLEIGH  DISTRIBUTION 


Y  =  -  14282. 


MOE 

CEP  1 

MEAN 

99.79 

STD  DEV 

11.56 

MSE 

133.74 

GED  MEAN 

99.12 

VARIANCE 

133.70 

COEF  VAR 

0.12 

SKEWNESS 

0.15 

KURTOSIS 

0.10 

RANGE 

70.75 

MINIMUM 

67.44 

.10  QUANTILE 

84.81 

.25 

92.04 

.50 

99.48 

.  75 

107.28 

.90 

114.56 

MAXIMUM 

138.19 

In  U  +  100 


/\  As. 


2 

CEP  3 

CEP  4 

79 

95.85 

108 

58 

12.64 

38.3 

09 

176.96 

1524 

11 

94 , 78 

100 

05 

159.74 

1460 

12 

0.13 

0.35 

15 

-  0.56 

0.36 

08 

4,59 

0.34 

25 

135.40 

248 

46 

7.73 

20.4 

78 

81.18 

59.0 

00 

87.60 

82.7 

48 

96.18 

105 

12 

103.97 

132 

59 

110.99 

156 

71 

143.13 

268 

6809 

/V 

CEP 

99 

11 

134 

99 

134 

0 

0 

0 

70 

67 

84 

92 

99 

107 

114 

137 
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CASE  V.  KEinULL  PISTKIIUJTION  w  •  ^ 

Y  -  (-  1442C050-J  .  I  In  U)  •* 


)E 

CEP  1 

CEP  2 

CEP  3 

CEP 

:an 

86.18 

85.80 

S-l  .25 

100 

’D  DEV 

5.23 

5.19 

11.47 

19. 

;e 

218.34 

228.61 

379.57 

382. 

:d  mean 

86.03 

85.64 

98. 

\RIAiNCE 

27.35 

26.97 

131.51 

382 

OEF  VAR 

0.06 

0.06 

0.14 

0. 

iCEU’NESS 

0.029 

0.028 

-  5.47 

-  0. 

j  IRTOSIS 

0.28 

0.27 

38.40 

0. 

vNGE 

32.67 

32.46 

112.79 

121 

'  NIMUM 

71.49 

71.14 

0.0 

420 

10  QUANTILE 

79.37 

79.04 

78.33 

74  . 

25 

82.87 

82.52 

81.89 

88  . 

50 

86.29 

85.89 

85.42 

100. 

75 

89.50 

90.12 

88.69 

113 

90 

92.50 

92.07 

91.56 

124 

1  .XI  MUM 

104.16 

103.61 

112.79 

163 

A 

6 

00 

1 

20 

24 

17 

2 

5 

0 
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CASE  VI.  WEIBULL  DISTRIBUTION  w  =  1.5 

2 

Y  =  (-  1442695.04  In  U) ^ 


MOE 

CEP  1 

CEP  2 

CEP  3 

CEP  4 

MEAN 

89.58 

89.28 

87.51 

102 

STD  DEV 

7.07 

7.04 

12.85 

25.7 

MSE 

158.65 

164.46 

321.29 

664 

GED  MEAN 

89.30 

89.00 

98.4 

VARIANCE 

50.07 

49.54 

165.09 

660 

COEF  VAR 

0.08 

0.08 

0.148 

-  0.25 

SKEVVNESS 

0.07 

0.07 

-  4.27 

-  0.04 

KURTOSIS 

0.22 

0.21 

27.91 

0.11 

RANGE 

44.55 

44.26 

125.57 

161 

MINIMUM 

69.36 

69.10 

0.0 

31.5 

.10  QUANTILE 

80.27 

79.99 

78.95 

67.6 

.25 

85.39 

85.12 

83.80 

86.4 

o 

LO 

• 

89.59 

89.25 

88.60 

101 

.75 

94.37 

94.07 

93.44 

120 

.90 

98.31 

98.00 

97.48 

133 

MAXIMUM 

113.91 

113.37 

125.57 

193 
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CASE  VII 


WEIBULL  DISTRIBUTION  w=  .5 


MOE 

CEP  1 

MEAN 

166.57 

STD  DEV 

42.61 

MSE 

6246.94 

GED  MEAN 

161.26 

VARIANCE 

1815.38 

COEF  VAR 

0.26 

SKEWNESS 

0.62 

KURTOSIS 

0.53 

RANGE 

250.79 

MINIMUM 

75.81 

.10  QUANTILE 

116.20 

.25 

135.54 

.50 

163.59 

.75 

191.58 

VO 

o 

220.85 

MAXIMUM 

326.60 

In  U)  ^ 


/\ 


2 

CEP  3 

CEP  4 

05 

87.45 

146 

99 

31.55 

99.2 

96 

1153.15 

11956 

32 

81.72 

113 

96 

995.65 

9840 

26 

0.36 

0.68 

70 

0.76 

1.39 

76 

1.48 

3.26 

53 

218.72 

746 

48 

18.67 

3.24 

74 

49.30 

41.2 

27 

66.56 

75.8 

51 

84.41 

126 

45 

105.80 

194 

60 

127.69 

282 

01 

237.40 

750 

4.27 

/s 

CEP 

170 

44 

6930 

164 

2023 

0 

0 

0 

258 

78 

117 

137 

165 

196 

227 

337 
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CASE  VIII.  V7EIBULL  DISTRIBUTION  w=  .1 

Y  =  (-  3.623886  In  U) 


MOE 

CEP 

1 

CEP  2 

CEP  3 

CEP  4 

MEAN 

2.6 

10  5 

3.6  10^ 

131.06 

5.1  lo"^ 

STD  DEV 

7.8 

10  5 

1.1  10^ 

335.02 

2.8  10^ 

MSE 

6.7 

lO^l 

1.4  10^^ 

1.13  10^ 

8.2  10^° 

GED  MEAN 

6.0 

lo"' 

7.4  lo"^ 

85.89 

674 

VARIANCE 

6.0 

io“ 

12 

1.3  10-^^ 

1.1  10^ 

7.9  10^° 

COEF  VAR 

2.97 

3.14 

2.56 

5.43 

SKEVJNESS 

7.95 

8.38 

14.16 

14.4 

KURTOSIS 

81.28 

89.17 

232.51 

255 

RANGE 

9.8 

10 

1.4  10^ 

6214.98 

5.4  10^ 

MINIMUM 

461. 

12 

476.59 

5.65 

7.2  10^ 

.10  QUANTILE 

6859 

.66 

7661.19 

30.28 

4.41 

in 

(N 

• 

1.8 

lo"^ 

2.2  lo"^ 

54.43 

77.0 

o 

in 

• 

5.9 

lo"^ 

7.4  lo'^ 

93.37 

1070 

.75 

1.7 

10^ 

2.2  10^ 

140.28 

9280 

.90 

5.6 

10^ 

7.3  10^ 

195.25 

8.3  lo"^ 

MAXIMUM 

9.8 

10^ 

1.4  10^ 

6220.63 

5.4  10^ 
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C0r4P ARISON  OF  ESTIMATORS 


USING  MEAN  SQUARED  ERROR  AS  A  CRITERION; 
MSE[CEP]  =  +  [E(CEP)  -  CEP]^ 


CASE  NO. 

DISTRIBUTION 

CEP  1 

A 

CEP  2 

A 

CEP  3 

\ 

CEP  4 

I 

RAYLEIGH 

136 

136 

165 

1549 

II 

UNIFORM 

109 

107 

368 

1744 

III 

POWER 

2460 

83,000 

444 

2279 

IV 

TRANSLATED 

RAY. 

134 

134 

177 

1524 

V 

WEIBULL  w 

2. 

218 

228 

379 

382 

VI 

WE I BULL  W 

1.5 

158 

164 

321 

664 

VII 

WEIBULL  w 

=: 

.5 

6246 

6930 

1153 

11956 

VIII 

V7EIBULL  w 

— 

.1 

6.77  10^^ 

6.78  l(f^ 

1.13  1(? 

8.17  10 

10 
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VII.  CONCLUSION 


In  general,  we  found  the  maximum  likelihood  power 
estimator  cSp  3  is  better  than  moment  method  power  estimator, 
in  that  it  does  not  have  the  limitations  and  exists  for  all 
parent  distributions  we  considered.  Using  minimum  mean 
squared  error  as  a  criterion,  the  estimators  based  on  the 
power  model  have  inferior  performance  in  the  "shallow  tailed" 
cases,  but  were  best  in  the  heavy  tailed  cases,  this  is 
especially  true  for  the  w  =  .1  Weibull  case,  in  which  it 
has  less  bias  than  its  competitors  and  has  relatively  small 
variance . 

We  found  the  Rayleigh  based  CEP  estimators  (CEP  1)  to 
be  possibly  the  best  overall,  especially  for  the  shallow 
tailed  cases.  The  maximum  likelihood  jackknife  Rayleigh 

/s 

estimator  (CEP  2)  had  no  better  performance  than  the  unbiased 
Rayleigh  estimator. 

The  value  of  the  sample  median  estimator  (CEP  4)  can 

be  found  for  small  samples  with  such  rapidity  that  in  certain 

cases  the  time  saved  may  compensate  for  the  accuracy  lost. 

In  practice,  when  analysts  find  several  impact  points 

comparatively  far  away  from  target,  they  should  choose 

/\ 

maximum  likelihood  power  estimator  (CEP  3) .  Otherwise  they 
could  choose  Rayleigh  based  CEP  estimator  (CEP  1) . 
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APPENDIX  1 


MAXIMUM  LIKELIHOOD  POWER  ESTIMATOR 


DIMENSION  Y( 200) , S ( 200 1 ) , A ( 200 1 ) 

DIMENSION  CEP(5G0) 

CALL  CVFLOW 

ISEED=470723 

READ  (5,100)  AO, SO 

READ  (5,200)  Ni. IST0P,S(1) ,A( 1) 

SN=FLOAT (MS) 

FA=-1,0/A0 

KS=500 

DC  1  KI=1,  KS 
C 

C  GENERATE  POWER  SQUARE  DEVIATE 

DC  5  i\=l,'NlS 
Y(N)=C.C 

CALL  RANDOM  (ISEED,U,1) 

Y(N)=SC'--(U^'=FA-l  .0  I 
5  CCMTNUE 
C 

C  SOLVE  NCNLINEAR  EQUATICN  UNDER  NEWTON'S  SECOND  METHOD 

DO  10  I=1,ISTOO 
ITERS=I 
A ( I+l ) -0. 0 
S( I+l)^C.O 

CALL  SF  (S(I  )  ,Y,NS,B) 

CALL  CF  ( S( T  ) , Y,NS  ,C  ) 

CALL  DF  (S(I ),Y,NS,0) 

E=alCG(  SID) 

UX  =  Si\/A  (  I  I+tr'^E-E 

VX  =  SN-A(I  )/S (I )-(A(I )  +  1.0)4D 

WX  =  -SN/  (  A(  T  )*:?2  ) 

xx  =  sr:/s  ( I  ) -0 

ZX  =  -SN-'A(  I  )/($(!  )T>?=2  )  +  (  M  I  )  +1.0)  YC 
0£T  =  WX’:^ZX-XX=>'=*2 
C 

C  CHECK  THE  DE Tr o M I NA T I  ON 

IF  (  GET  .cQ.' .'■•  )  CO  TO 
UU=ZX*UX-XX-VX 
VV=WX*VX-XX^UX 

c 

C  WHETFiER  TO  STOP  "HE  ITERATION 

CHI=-(  LX'AUUi-VX-VV  ) /DET 
I F  (CHI .LE.0.C4U )  GO  TO  03 
A (  I  +  l )  =  A( I  )-UU/OFT 
S(  I  +  l  )  =  S( I  )-VV/DcT 
IF  (S(I+1  )  .GT.G.O)  GO  TO  10 
S (  I  +  l  )  =  S(  I  )/2.0 
10  CCNTINLIE 
35  ACUT=A( ITERS ) 

.SGUV=S(  ITERS  ) 

C 

C  ESTIMATE  CEP 

AK  =  2-*(  1.  O/AOUT  )  -1.0 

CEP(KI  )  =  (SOUT=!=AK  )*=?C..5 

WRITE  (6,40J)  SOUT, AGUT, ITERS, C5P(KI) 

1  CONTINUE 

CALL  HISTG  (CEP,KS,0) 

99  STOP 
END 


SUEROUTINE  PF  (S,Y,NS,S) 
DIMENSION  Y(200) 

B=C.O 

DC  30  N=1,NS 
B=e+ALOG( S+Y(N ) J 
30  CCN'TINUE 
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RETURN 

END 

SU6RCUTIN5  CF  (5,Y,NS,C) 
DIMENSION  Y(200) 

C  =  0.0 

DO  AO  N-1,NS 
C=C+1. 0/(5+Y(N) 

AC  CONTINUE 
RETURN 
END 

SUEROUTINE  DF  (S,Y,NS,D) 
DiNENSICN  Y(200) 

0  =  0.0 

no  50  N  =  1,.NS 
D=D+1.0/( S+Y(N ) ) 

50  CCNVINUE 
RETURN 
END 
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APPENDIX  2:  RAYLEIGH  UNPIASEO  AND  JACKKNIFE  ESTIMATOR 


DIMENSION  Y(50,50),SP<50),PSE{50), CEP{ 500) 

DIMENSION  RAY (500) 

C 

C  GENERATE  RAYLEIGH  DEVIATE 

CALL  OVFLOW 
ISEED=470722 
C 

C  READ  IN  ORIGINAL  POP.  PARAMETER  VC  ,  TOTAL  SAMPLE  SIZE 

C  NS,  SUBGROUP  SAMPL^  SIZE  KS 

READ  (5,100)  VJ,NS,KS 
CPO^l.  17741C02'':‘S0RT(  VO  ) 

MS=NS/KS 
SM=FLOAT{MS ) 

SN=FL3AT(NS) 

SK=FLGAT{KS) 

WRITE  (o,2C0)  VD, C?0,NS,MS ,KS 
KE=5C0 

DO  1  KCE-1,KF 
DO  5  m=i,m;s 
DO  5  K=1,KS 
Y(M,K)=0,0 

CALL  RANDOM  ( ISEED,U, 1 ) 

Y  (  M  ,  K  )  =  iOOOO  .  0=''  {  1  .  U/U-  1 . 0  ) 

5  CONTINUE 
C 

C  ESTIMATE  THE  PARAMETER  BASED  ON  ALL  DATA 

T  =  C.O 

DO  10  M  =  1,MS 
DC  IG  K=1,KS 
T  =  T+Y(yi,K) 

10  CONTINUE 

TM  =  SCRT(T/(SN'-'=2.0)  ) 

C 

C  MAXIMUM  LIKELIHOOD  UNBIASED  ESTIMATER 

RAY(KKE  )=1.184S*TM 
C 

C  ESTIMATE  THE  PARAMETER  BASED  ON  ALL  DATA.  EXCEPT  THOSE 

C  IN  THE  ITH  GROUP 

DO  20  M=1,MS 
S  =  0.0 

DC  25  1=1, MS 
IF  (I .EQ.M)  G3  TO  25 
DC  26  k=1,KS 
S=S+Y( I ,K) 

26  CCN'TINJE 
25  CONTINUE 

SP(M)=SQRT(S/( (SN-SK)*2.0) ) 

20  CONTINUE 
C 

C  COMPUTE  THE  MS  PSEUDO  VALUES 

DC  4C  M  = 

PSE(  M)  =  SMMTS-(  SM-1.0)=;=SP(M) 

40  CONTINUE 
C 

c  estimate  jackknife  parameter 

ACK=C.O 
DO  45  M=1,MS 
ACK=ACK+PSE(M) 

45  CONTINUE 

VJACK=ACK/SM 
S  C  K = 0  .  o 

c 

C  COMPUTER  TfHE  CONFIDENCE  INTERVAL 

00  50  1=1, MS 

SCK=SCK+( PSE( I )-VJACK)**2 
50  CONTINUE 
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SyJA=SCK/  (  SM=I‘{  ?M-1 ,0)  ) 

UPE  =  VJi  CK-2.  D-SOra  (SVJM 
PLB  =  VJACK  +  2.0=?SQKT  (SVJ4) 

C 

C  ESTIMATr  C5?  FR3M  j:.CKKNIFS  PARAMETER 

CEP(KRF  I=1.1774*VJAC!< 

WRITE  (6,303)  V J AC K , UP B , PL  6, CE P ( Kk E ) 

1  CCNTIiNUE 

CALL  HISTG  (RAY,KE,C) 
call  HISTG  (CFP,i<r,G) 

STOP 

END 
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